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Abstract. We report on an ongoing investigation into a seismic calibration of solar models 
designed for estimating the main-sequence age and a measure of the chemical abundances of 
the Sun. Only modes of low degree are employed, so that with appropriate modification the 
procedure could be applied to other stars. We have found that, as has been anticipated, a 
p^ ' separation of the contributions to the seismic frequencies arising from the relatively smooth, 
glitch-free, background structure of the star and from glitches produced by helium ionization 
and the abrupt gradient change at the base of the convection zone renders the procedure more 
robust than earlier calibrations that fitted only raw frequencies to glitch-free asymptotics. As 
fsj . in the past, we use asymptotic analysis to design seismic signatures that are, to the best of 
rsj ' our ability, contaminated as little as possible by those uncertain properties of the star that 
are not directly associated with age and chemical composition. The calibration itself, however, 
employs only numerically computed eigenfrequencies. It is based on a linear perturbation from 
a reference model. Two reference models have been used, one somewhat younger, the other 
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somewhat older than the Sun. The two calibrations, which use BiSON data, are more-or-less 
consistent, and yield a main-sequence age t@ — 4.68 ± 0.02 Gy, coupled with a formal initial 
heavy-element abundance Z — 0.0169 ± 0.0005. The error analysis has not yet been completed, 
C/3 , so the estimated precision must be taken with a pinch of salt. 

Keywords. Sun: helioseismology, Sun: abundances, stars: abundances, stars: oscillations, stars: 
fundamental parameters. 
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1. Introduction 



(*C) • The only way by which the age of the Sun can be estimated to a useful degree 



of precision is by accepting the basic tenets of solar-evolution theory and measuring 
those aspects of the structure of the Sun that are predicted by the theory to be in- 
dicators of age. The structure measurements must be carried out seismologically, and 
evidently one expects greatest reliability of the results when all the available helioseis- 
^- , mic data are employed. However, the most relevant modes are those of lowest degree, 
JJ~j ■ because it is they that penetrate most deeply into the energy-generating core where 
the relic helium-abundance variation records the integrated history of nuclear trans- 
mutation. Moreover, it is also only they that can be measured in other stars. There- 
fore, there has been some interest in calibrating theoretical stellar models using only 
low-degree modes. The prospect was first discussed in detail by Christensen-Dalsgaard 



(1984, 88), Ulrich (1986) and Gough (1987) although prior to that it had already been 



pointed out that the helioseismic frequency data that were available at the time indi- 
cated that either the initial helium abundance Yq or the age £©, or both, are some- 
what greater than the generally accepted values (Gough 1983), an inference which is 
consistent with our present findings. Subsequent, more careful, calibrations were car- 
ried out by |Guenther (1989)| |Gough fc Novotny (1990)} |Guenther k Demarque (1997) 
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Weiss fc Schlattl (1998) Dziembowski et al. (1999) Gough (2001) and Bonanno, Schlattl 



& Paterno (2002). Not all of them addressed the influence of uncertainties in Yq on the 
determination of i Q . 

As a main-sequence star ages, helium is produced in the core, increasing the mean 
molecular mass /i, preferentially at the centre, and thereby reducing the sound speed. 
The resulting functional form of the sound speed c(r) depends not only on age i Q but 
also on the relative augmentation of //(r), which itself depends on the initial absolute 



value of /x, and hence on Yq. Gough (2001) tried to separate these two effects using the 



degree dependence of the small separation d n ,i — 3(2^ + 3)~ l (v n .i — Vn-1,1+2) of cyclic 
frequencies v n ,i, where n is order and I is degree. This is possible, in principle, because 
modes of different degree and similar frequency sample the core differently. However, that 
difference is subtle, and the sensitivity to the relatively fine distinction between the effects 
of £q and Yq on the functional form of c(r) in the core is low. Consequently the error in 
the calibration produced by errors in the observed frequency data is uncomfortably high. 
This lack of sensitivity can be overcome by using, in addition to core-sensitive seismic 
signatures, the relatively small oscillatory component of the eigenfrequencies induced by 
the sound-speed glitch associated with helium ionization (Gough 2002), whose amplitude 
is close to being proportional to helium abundance Y (Houdek & Gough 2007a). The 
neglect of that component in the previously employed asymptotic signature d nt i had 
not only omitted an important diagnostic of Y but had also imprinted an oscillatory 
contamination in the calibration as the limits (fci, fe), where k = n + ^l, of the adopted 
mode range was varied (Gough 2001). It therefore behoves us to decontaminate the core 
signature from glitch contributions produced in the outer layers of the star (from both 
helium ionization and the abrupt variation at the base of the convection zone, and also 
from hydrogen ionization and the superadiabatic convective boundary layer immediately 
beneath the photosphere) . To this end a helioseismic glitch signature has been developed 



by Houdek & Gough (2007a) , from which frequency contributions 5v r ^i can be computed 
and subtracted from the raw frequencies v n j to produce effective glitch-free frequencies 
v SUt i to which a glitch-free asymptotic formula (|2.10|) can be fitted. The solar calibration 
is then accomplished by interpolating the theoretical seismic signatures computed on a 
grid of solar models to the observations, using a standard grid to compute derivatives 
with respect to i Q and Yq, and a carefully computed reference solar model designed to 
be close to the Sun. The result of the first preliminary calibration by this method, using 



BiSON data, has been reported by Houdek fc Gough (2007b)| Here we enlarge on our 



discussion of the analysis, and we augment our results with a calibration based on a 
second reference solar model. 

2. The seismic diagnostic and calibration method 

Any abrupt variation in the stratification of a star (relative to the scale of the inverse 
radial wavenumber of a seismic mode of oscillation), which here we call an acoustic 
glitch, induces an oscillatory component in the spacing of the cyclic eigenfrequencies of 
seismic modes. Our interest is principally in the glitch caused by the depression in the 
first adiabatic exponent 71 = (<91np/<91np) s (where p, p and s are pressure, density and 
specific entropy) caused by helium ionization, which imparts a glitch in the sound speed 

c{r). The deviation 

5vi := v t - v si , (2.1) 

where i := (n,l), of the eigenfrequency i/j from the corresponding frequency i/ S i of a 
similar smoothly stratified star is the indicator of Y that we use in conjunction with the 
indicators of core structure to determine the main-sequence age. 

Approximate expressions for the frequency contributions Svi arising from acoustic 



Progress report on solar age calibration 











; Sun (BiSON) 






Model S 


A A* njui il 






A «A*V/\»V ♦ 










M I V " 






, fly ; 


| J 
























\* * AM 






u 


* * 








«* % 4 ? 






* 


r * A * * 






** ? 


' * , 






5 * * * 


r x » „ * 






* * 










•a 






yi *„**«** 


r ^^ 






*** 


♦ 






*J£» 
















♦* 








♦ J 
















* ; 


: ..... 






. 



2000 3500 3000 3500 4000 1000 1500 2000 3500 3000 3500 



' (/"Hz) 



' (/*Hz) 



Figure 1. Top left: The symbols (with error bars obtained under the assumption that the raw 
frequency errors are independent) represent second differences, A2i^, of low-degree solar fre- 
quencies from BiSON. Top right: The symbols are second differences A2^ o f adiabatic pulsation 
eigenfrequencies of solar Model S of |Christensen-Dalsgaard et al. (1996)| The solid curves in 
both panels are the diagnostics (12. 2[) - (|2.8[) . whose eleven parameters have been adjusted to fit 
the data optimally. Bottom: The symbols denote contributions 5v to the frequencies produced 
by the acoustic glitches of the Sun (left panel) and Model S (right panel) . 



glitches in solar-type stars were recently presented by Houdek & Gough (2007a) Here 



we improve them by adopting the appropriate Airy functions Ai(— x) that are used as 
comparison functions in the JWKB approximations to the oscillation cigenfunctions, as 



did Houdek & Gough (2007b) The complete expression for 8vi is then given by 



5vi = 8-yVi + 8, 



7^2 



(2.2) 



where 
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arises from the variation in ji induced by helium ionization, and 
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(2.4) 



arises from the acoustic glitch at the base of the convection zone resulting from a near 
discontinuity (a true discontinuity in theoretical models using local mixing-length theory 
with a non-zero mixing length at the lower boundary of the convection zone) in the second 
derivative of density. Here, m = 3.5 is a constant, being a representative polytropic index 
in the expression for the approximate effective phase ip appearing in the argument of the 
Airy function, and /?, fj and jl are constants of order unity which account for the relation 
between the acoustic glitches caused by the first and second stages of ionization of helium 
(Houdek & Gough 2007a); r is acoustic depth beneath the seismic surface of the star, and 
T ~ \J1vq is the total acoustic radius of the star; An and t\\ are respectively the acoustic 
width of the glitch and its acoustic depth beneath the seismic surface. The argument of 
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the Airy function is x = sgn(-0)| 3-0/2 | 2 / 3 , where 

ip(r) — kujt — (m + 1) cos - [(to + l)/u>f] if f > r t , (2-5) 

and 

ip(r) = \k\ujt- (to + 1) ln[(m + l)/wf+ |«|] iff<r t , (2.6) 

in which f = r + u; en, with oj = 2irv, and rt is the location of the upper turning point 

) 2 /oj 2 f 2 ] 1 / 2 , and Ki = K(rJTi_i) and kh — k(tli). In 

1) tos -1 [(to + l)/f c w] + tt/4, (2.7) 

ii, 7ii) eii) A c , r c , e c ), a = 1, ...,7, are found by 

3 

A 2i ^ = i/„_i,j - 2z/ n ,i + t/ n+ i,i ~ A 2i (5 7 ^ + <5 C ^) + 2J afc^r* — 9i( v ]^Va) (2.8) 

fc=0 

to the corresponding observations by minimizing 

£ g - (A 2i */ - gJC^fajV - 9j ) (2.9) 

using the value of z^o obtained by fitting to (|2.10p . where C^J- is the (i,j) element of 
the inverse of the covariance matrix Ca of the observational errors in A 2 ^, computed, 
perforce, under the assumption that the errors in the frequency data v>i are indepen- 
dent. The last term in equation (|2.8p approximates smooth contributions arising, in part, 
from wave refraction in the stellar core, from hydrogen ionization and from the supera- 
diabaticity of the upper boundary layer of the convection zone, introducing four more 
fitting coefficients cik = i] a , k = 0, ..., 3, a = 8, ..., 11. The covariance matrix C^ Q7 of the 
errors in r\ a were established by Monte Carlo simulation. 

The outcome of the fitting to the BiSON data (Basu et al. 2007) and to the adiabati- 
cally computed eigenfrequencies of solar Model S (Christensen-Dalsgaard et al. 1996) is 
displayed in Figure 1: the upper panels display the second differences, together with the 
fitted formula (|2.8|) . the lower panels display the corresponding contributions 5vi to the 
frequencies of oscillation from the acoustic glitches. All the frequencies displayed in the 
figure have been used in equation (|2.9| for fitting l|2.8|l . 

To the resulting glitch- free frequencies v S i, derived from equation (|2.ip , of both the 

solar observations and the eigenfrequencies of the reference solar model, is fitted the 

asymptotic expression 

. . „ AL 2 -B 2 CL 4 -DL 2 +E 4 FL 6 -GL 4 + HL 2 -I 6 . . . 
Vsi^in+^l+ejVQ v 3 u g v a = s(v si ;t,p), 

(2-10) 

by minimizing {y S i — Si)C si ^ (vj — Sj), where L 2 = 1(1 + 1) and C s is the covariance 
matrix of the observational errors in i/ si , from which we obtain both the coefficients 
£/3 = (^o, £, A, B, C,D,E, F, G, H, I), (3 = 1, ..., 11, and the covariance matrix Ccgs of 



the errors. Following Gough (2001) we carry out this fitting in the frequency range given 
by k\ ^ k ^ &2, where k = n + ^l and ^ / ^ 3, and we vary k\ and fc 2 . Each of the 
parameters ^g represents an integral of a function of the equilibrium stratification. The 
integrals A, C and F are of particular importance to our analysis, because C and F are 
dominated by conditions in the core, and, although the contributions to A from the core 
and the rest of the star are roughly equal in magnitude (and potentially have opposite 
signs), the latter is relatively insensitive to £© and Y$. The integrands in the remaining 
integrals are either more evenly distributed throughout the Sun or are concentrated near 
the surface. 
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Table 1. Partial derivatives H a j obtained from two sets of calibrated evolutionary models for 
the Sun. Values with respect to age £© are in units of Gy _1 . 

(dA/dt ) z (dA/dZ) tQ (dC/dt Q )z (dC/dZ) tQ [0(-<J7i/7i)/0t Q ]z [d(-S-yihi)/dZ\ te) 
-0.0469 -0.584 0.677 36.8 -0.00656 0.442 

We have carried out age calibrations using combinations of the parameters 

C a = (i4,C,-*7i/7i), a = 1,2, 3, (2.11) 

where —571/71 = An/ \/Q/kvqAii is a measure of the maximum depression in 71 in the 
second helium ionization zone. Presuming, as is normal, that the reference model is 
parametrically close to the Sun, we consider the reference value Q to be approximated 
by a two-term Taylor expansion of £ a about the value £® of the Sun: 



«^»-l£J z Afe -lafJ, e AZ + E <»' (2 ' 12 > 

where A i© and AZ are the deviations of age t© and initial heavy-clement abundance 
Z from the reference model, and e^ a are the formal errors in the calibration parameters 
whose covariance matrix C^ a p can be derived from C^ps and C r?Q7 . A (parametrically 
local) maximum-likelihood fit then leads to the following set of linear equations: 

H a jC ( ~ al3 Hpk < dok — H a jC^ al3 A j3 , (2-13) 

in which Ofc = (At©, AZ) + eefe = Oofe + eefe, k = 1,2, is the solution vector subject 

H aj = [(d( a /dt Q ) z , (dU9Z) t@ ]J = lf2. 

A similar set of equations is obtained for the formal errors eQk' 

H ajC^ a/3 Hp k e ek = H aj C^ af3 e C p , (2.14) 

from which the error covariance matrix Cskq = eefcCeg can be computed from C^ a p. 
The partial derivatives H a j were obtained from the two sets of five calibrated evolution- 



to (correlated) errors eefc, A^ — C,p — Cp + e c/3 = ^0/3 + e (p, and the partial derivatives 



ary models for the Sun that were used in a similar calibration by Houdek & Gough (2007b) 
computed with the evolutionary programme by Christensen-Dalsgaard (1982)| and adopt 



ing the Livermoore equation of state and the OPAL92 opacities. One set of models has 
a constant value for the heavy-element abundance Z = 0.02 but varying age; the other 
has constant age but varying Z. Note that, for prescribed relative abundances of heavy 
elements, the condition that the luminosity and radius of the Sun agree with observation 
defines a functional relation between Yq, Z and £©. The values of the partial derivatives 
H a j are listed in Table 1 . 

3. Results 

To illustrate the effect of taking 5vi into account, we compare in Figure 2 a first as- 
sessment of A(k\,k2) using the glitch-frce frequencies v S i (left panel) with that obtained 
from the raw frequencies Vi (right panel). Recall that A represents a functional of the 
equilibrium structure of the star, and should not vary with k\ and k^. The range of values 
for A is the lower for v S i, as we had anticipated. We believe that the upturn of A for 
low values of k\ in the left panel of the figure is a result of the failure of the asymptotic 
formula (|2.2[) - (|2.4p when v h is low. The dipping of A at high values of k\ and low values 
of ki occurs because the frequency range is too small for a reliable determination of the 
fitting coefficients ^g. We therefore adopt intermediate values for k\ and high values for 
/c2, for which A is insensitive to the selected frequency range. 
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Figure 2. Asymptotic fitting coefficient A (see equation I2.10|l as functions of fci and fe 
(k — n + i-Z). Results are shown for fitting (|2.1Up to the glitch-free frequencies v si (left panel) 
and to the raw frequencies Vi (right panel). 



Age calibrations using different combinations of the parameters £a and two different 
reference models are summarized in Table 2. The younger reference model is 'Model S' 
(Christensen-Dalsgaard et al. 1996) which has age i Q = 4.6 Gy; the second is 'Model T", 
which has age i© = 4.7 Gy. The same physics was adopted in the evolutionary calculations 
of both models. We notice in Table 2 that the calibration for the combination (A, C), i.e. 
without £71/71, is less stable to a change in the reference model than are the calibrations 
with combinations in which £71/71 is included, and therefore is less reliable, as we have 
explained in the introduction. If we ignore in Table 2 the results for {A, C) and combine 
the others, we obtain 



t. 



4.68 ± 0.02 Gy, Z = 0.0169 ± 0.0005. 



Including the calibrations with (A, C) does not change the outcome. Error contours 
corresponding to the calibration from Model S in the first row of Table 2 are plotted in 
Figure 3. Corresponding contours for Model T arc the same, except that their centres are 
displaced to (4.677 Gy, 0.0170). One can adduce from our description of the analysis in 
Section 2 that our current treatment of the errors is not completely unbiassed; however, 
the potential bias is of the order of only \6vi/vi\, which is small. 

The age we have found is greater than currently accepted values. The values of Z 
are somewhat smaller than those of Models S and T (0.01963), but we hasten to point 
out that they should not be regarded strictly as statements about the initial heavy- 
element abundance, but rather as measures of the opacity in the radiative interior. 



Asplund et al. (2004) have argued that the photospheric abundances of C, N and O 



had previously been overestimated, suggesting that the actual total heavy-element abun- 
dance is even lower than had previously been believed. However, that cannot imply that 
the opacity in the solar interior is necessarily comparably lower because it has been im- 
plicitly calibrated here (by accepting the tenets of solar-evolution theory, and the OPAL 
opacity calculations upon which the models are based), and indeed the opacity has al- 
ready been determined seismologically from a broader spectrum of modes than has been 
adopted here (Gough 2004). The matter raised by Asplund etal. therefore challenges 
either the opacity calculations, the nuclear reaction rates, or the basic physics of stellar 
evolution, not hclioseismology as some spectators have surmised. As we know already 
from seismological structure inversions, the solar models are not accurate by helioseis- 
mological standards. Therefore the properties inferred from these calibrations could be 
more contaminated by systematic error than by errors in the observed frequencies. 
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Table 2. Age calibrations with different combinations of £ a and for the two reference models: 
ModelS with an age t© = 4.6 Gy and Model T with an age £q = 4.7 Gy. The first two columns 
show the results adopting Model S as the reference model, the third and fourth columns display 
the results for Model T. 
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Figure 3. Error ellipses for the calibration 
using all three parameters £ a and Models S as 
the reference model: solutions (£q , Z) satisfy- 
ing the frequency data within 1, 2 and 3 stan- 
dard errors in those data reside in the inner, 
intermediate and outer ellipses, respectively. 
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Discussion 

Christensen-Dalsgaard: A comment: with SONG we expect to be able to carry out 
a similar analysis of distant stars, on which we of course know much less a priori. 

Houdek: This seismic diagnostic has been developed with the aim to be able to use it 
also for distant stars. The accuracy of the observed frequency data required for such a 
diagnostic analysis is one part in 10 4 or better. 

S. Vauclair: Two small comments which are actually more relevant for stars that are 
slightly more massive than the Sun; First, I would like to point out that in case of 
helium settling below the convective zone the effect of the helium gradient in the second 
differences may become more important than the convective border, and than the effect 
of helium ionization. Second, the so-called asymptotic theory, which is very useful, may 
become quite wrong in some cases, at the end of the main sequence or the beginning 
of the sub-giant branch. The small frequency separation, which is always positive in the 
asymptotic theory, can become negative. 

GOUGH & Houdek: You are certainly correct in implying that the amplitude of the 
oscillatory contribution to the second differences arising from helium settling beneath the 
convection zone can be greater in stars more massive than the Sun, which have shallower 
convection zones, although whether or not it is more important than the ionization 
signature depends upon the issue in question. The cumulative amount of settling increases 
with time, and therefore is a potential indicator of age. But the sound-speed profile that 
it produces depends on uncertain fluid-dynamical issues associated with the tachocline, 
the recession of the convection zone, and possible overshooting, so we would be wary of 
attempting to use its seismic signature in an age calibration. In the current state of our 
understanding we would instead prefer to separate it from the ionization signature and 
then ignore it, as we have for the Sun; that course is possible provided that the helium 
ionization zone is acoustically far from the base of the convection zone. We would use it 
separately to investigate tachocline structure, however, as indeed we are in the process 
of doing for the Sun. 

One cannot deny that conditions in some stars might be such as to render it impossible 
to develop an adequate asymptotic theory of low-degree acoustic modes, although we do 
not share your apparently implied pessimism. It is perhaps worth pointing out that there 
have been instances when an asymptotic formula developed in one set of circumstances 
has been misused by applying it without modification in another, in which the conditions 
for the validity of the theory are not satisfied; if that is what you mean by "so-called 
asymptotic theory" we must surely agree. We must point out, however, that it is not true 
that even the simple asymptotic glitch-frce formula (|2.10[) precludes a negative so-called 
small frequency separation. 



